Backward correlations and dynamic heterogeneities: a computer study of ion 

dynamics 
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We analyse the correlated back and forth dynamics and dynamic heterogeneities, i.e. the presence 
of fast and slow ions, for a lithium metasilicate system via computer simulations. For this purpose we 
define, in analogy to previous work in the field of glass transition, appropriate three-time correlation 
functions. They contain information about the dynamics during two successive time intervals. First 
we apply them to simple model systems in order to clarify their information content. Afterwards 
we use this formalism to analyse the lithium trajectories. A strong back-dragging effect is observed, 
which also fulfills the time-temperature superposition principle. Furthermore, it turns out that 
the back-dragging effect is long-ranged and exceeds the nearest neighbor position. In contrast, the 
strength of the dynamic heterogeneities does not fulfill the time-temperature superposition principle. 
The lower the temperature, the stronger the mobility difference between fast and slow ions. The 
results are then compared with the simple model systems considered here as well as with some 
lattice models of ion dynamics. 

PACS numbers: 66.30.Dn 



I. INTRODUCTION 

One basic characteristics of amorphous ion conductors 
is the strong frequency dependence of the conductivity 
at sufficiently low temperatures ^, ||, 0] . According to 
hnear response theory, the dispersion in alv) is equivalent 
to a non-difiusive mean square displacement {r'^{t)) and 
thus to the presence of correlated back and forth jumps 
0. Neglecting possible dynamic correlations among ad- 
jacent ions the relation between both quantities reads 



dt{d/dt)w{t)exp{-i2TTiyt) (1) 



where q denotes the charge and p the density of the mo- 
bile ions. The function w{t) is defined as 



w{t) = {d/dt){r^{t)). 



(2) 



For very high frequencies one observes local dynamics 
which is strongly system-dependent. For example, in a 
sodium silicate system ^ one observes for v < 10^^ Hz 
a continuous decrease of cr(i^) with decreasing v which 
can typically be attributed to non-local dynamics |^ 
and becomes stronger for decreasing temperature. Ac- 
tually, this conclusion has been explicitly verified in re- 
cent simulations for lithium metasilicate . In a double- 
logarithmic representation at low temperatures the ap- 
parent exponent decreases from one to zero until the d.c. 
plateau is reached, i.e. (j{u) = Ud.c- 
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In general, the complexity of ion dynamics in disor- 
dered systems is, on the one hand, related to the static 
disorder of the material and, on the other hand, to the 
Coulomb interaction among the ions giving rise to dy- 
namic disorder. Therefore, one might expect that both 
the static as well as the dynamic disorder enhances the 
number of correlated back and forth jumps (see, how- 
ever, Ref. [^) although in recent simulations of sodium 
silicate at T = 2000 K no such correlations have been 
observed Furthermore, the presence of different envi- 
ronments of the network might imply that at a given 
time some ions are more mobile than other ions, i.e. 
there exist dynamic heterogeneities. Unfortunately, no 
direct information about dynamic heterogeneities can be 
gained from conductivity experiments or, equivalently, 
from the time dependence of the mean square displace- 
ment. Some progress has been achieved on the basis of 
simulations. Comparing the mean square displacement 
of different ions during a given time a broad distribu- 
tion has been observed [^. This basically corresponds 
to the presence of a large non-gaussian parameter which 
is known to represent dynamic heterogeneities quite well 
Other groups found an interesting spatial structure 
of the mobile regions |l^ . 

Another system with intrinsic complex dynamics is a 
glass-forming liquid. A tagged particle in a glass-former 
experiences a dramatic slowing down with decreasing 
temperature |Q. To a large extent this is related to 
the cage effect since at low temperatures the neighbour 
particles have a strong confining effect on the central par- 
ticle jl^ . This gives rise to strongly subdiffusive dynam- 
ics, i.e. to correlated back and forth dynamics. Fur- 
thermore glass-forming systems display dynamic hetero- 
geneities ^ ^ . This property has been quanti- 
fied by invoking appropriate three-time correlation func- 
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tions, containing information about the dynamics during 
two subsequent time intervals |2^, |2^, . 

Using these three-time correlation functions one can 
hope to answer several basic questions about the com- 
plexity of ion dynamics. How relevant are backward 
correlations? Are back and forth correlations restricted 
to nearest-neighbor ionic positions? How does the ten- 
dency of these back and forth correlations depend on the 
time scale of investigation? Do dynamic heterogeneities 
depend on temperature and thus invalidate the time- 
temperature superposition principle, observed for many 
other quantities like the conductivity? This and other 
aspects of ionic dynamics will be analysed in this paper 
by invoking appropriate three-time correlation functions. 
They will be applied to computer generated trajectories 
of lithium metasilicate (Li20)(Si02). Very recently, first 
measurements of three-time correlations have been con- 
ducted via multidimensional NMR. These results clearly 
showed that dynamic heterogeneities are omnipresent in 
disordered ionic conductors Q 

The organization of this paper is as follows. In Sec- 
tion 2 we describe the technical aspects of the simulation 
and discuss the numerical tools. Section 3 introduces the 
concept of three-time correlations. Section 4 contains 
the analysis of simple model systems for which we clar- 
ify the information content of the three-time correlation 
functions. The results of our simulations as well as their 
interpretation are presented in Section 5. We close with 
a discussion and a summary in Section 6. 




Figure 1: Sketch of a single-particle dynamics in order to 
clarify the definition of roi and ru- 

the lithium subsystem has been equilibrated before start- 
ing the production run. In particular at the lower tem- 
peratures the network fluctuations are very small (mean 
square displacement of 0.3 A^). 



III. THREE-TIME CORRELATIONS 
A. Definition of three-time correlations 



II. SIMULATION 

The potential energy for the lithium silicate system is 
chosen to be the sum of a Buckingham and a Coulomb 
pair potential denote the species lithium, oxygen, or 
silicon, respectively) 

U.^i^r) = eM-B.,r). (3) 

The parameters have been determined by Habasaki et 
al. JUJ, |2^. Details of our simulation can be found in 
10, |26| . Summarizing, we performed molecular dynam- 
ics simulations with a time step of 2 fs and a density 
of p — 2.34 g cm^"^. Periodic boundary conditions were 
used. The system size is 1152 particles, thus containing 
384 lithium ions. The trajectories were generated by an 
appropriately modified version of the MOLDY software 
package, supplied by K. Refson |27|. The length of the 
production runs was 16 ns after an equilibration time of 
ca. 10 ns at the two lowest temperatures. Even at T = 
640 K the mean square displacement during the produc- 
tion time was larger than 60 (nearest-neighbour Li- 
Li distance 2.6 A); see Ref. Q]. All configurations were 
equilibrated at T = 1500 K. The computer glass tran- 
sition is approximately 1100 K |2^. In this work we 
present simulations for five temperatures (T = 1240 K, T 
= 980 K, 750 K, 700 K, and 640 K). For all temperatures, 



As mentioned above, the dispersive behavior of the 
conductivity or, equivalently, the non-diffusive behavior 
of the mean square displacement, has been related to 
correlated back and forth dynamics. Here we want to 
introduce a formalism which allows one to elucidate the 
properties of back and forth dynamics in detail. We re- 
mind the reader that standard observables for the char- 
acterization of dynamical properties like the mean square 
displacement correlate the position of individual ions at 
two successive times to and ti. For a calculation of {r'^{t)) 
one averages over the configurations at all times to and ti 
such that ti—to = t. Back and forth dynamics, however, 
is related to the properties of particles during two subse- 
quent time intervals. Formally, this can be described as 
a three-time correlation, including a third time t2 > ti. 
The idea is sketched in Fig. |l]. roi denotes the distance 
a particle moves during the first time interval of length 
^01 = ti—to- The value of ri2 denotes the motion during 
the second time interval of length ti2 = t2 — ti as pro- 
jected on the direction of the motion during the first time 
interval. In case of a backjump, as shown for the exam- 
ple in Fig. |l|, the value of ri2 is counted negative. The 
additional information about the dynamics during two 
successive time intervals as compared to the dynamics 
during a single time interval is contained in the condi- 
tional probability function p(ri2|roi) which denotes the 
probability for a specific value of ri2 under the condition 
that the particle has moved the distance roi in the first 
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time interval. 



B. Moments 

Rather than analysing the full probability function, we 
concentrate on the first moment r(roi) and the second 
moment w(roi) = ((fi2 — r{rQi))'^{roi)). The interpreta- 
tion of both functions is straightforward, f(roi) contains 
information about the relevance of back and forth dy- 
namics. In case that the direction of the dynamics during 
two successive time intervals is uncorrelated one expects 
^(''oi) = 0. In contrast, a negative value of f is direct 
evidence of the presence of back and forth dynamics. 

The second moment w(roi) yields information about 
the presence of dynamic heterogeneities. In case that all 
particles have the same mobility, the distance moved in 
the second time interval is independent of the distance 
moved in the first time interval. Thus v(rQi) would not 
depend on roi. A dependence on rpi will be observed if 
there exist fast and slow ions. The subensemblc of ions 
with small roi will preferably contain slow ions whereas 
for ions with large roi it is vice versa. Therefore, ions 
with small rpi will on average move less in the second 
time interval ti2 than ions with large rpi, resulting in a 
monotonous increase of v{roi) with tqi. 



First moment in the limit toi 







It is possible to establish a direct relation between the 
mean square displacement and the first moment f(roi) for 
the case of a discrete hopping model in the limit ^oi 
0. For reasons of simplicity we take a ID model with 
distances do between the individual sites. One may start 
with the simple relation 

{r\tQi+h2)) = (r2(ioi)>+2(r(toi)r(ii2)) + (r2(<i2)) (4) 

which is valid for stationary processes. On the left side 
one can perform a linear expansion around ioi = 0. On 
the right side one may use the fact that for very small 
toi the term (r^(toi)) can be written as re//c?o^oi (cor- 
responding to the short-time diffusion in pure hopping 
models with an effective escape rate Fe//). Since for 
very short times the system can only jump to the nearest 
neighbor site the term {r{toi)r(ti2)) can be expressed as 
Teffdotoifido). Inserting these relations into Eq.^ one 
finally ends up with 



rjdo) 
do 



w{0) 



(5) 



This relation directly shows that exactly in case of dif- 
fusive dynamics, i.e. w{t) = const, one has f{do) — 0, 
i.e. no correlated back and forth dynamics. For subdiffu- 
sive behavior one obtains (beyond a possible oscillatory 
regime of the mean square displacement) < w{t) < 



w{0) and thus -1/2 < f(do)/do < 0. Thus subdiffu- 
sive behavior is equivalent to the presence of correlated 
forward and backward jumps. Furthermore validity of 
Eq.^ implies that the first moment f(do) has a lower limit 
— (l/2)(io which is reached if w{t) — 0, i.e. for the long- 
time limit of localized dynamics where any forth jump 
is followed by a back jump . Finally one can see that 
the range of dispersion, i.e. w{t — s- oo)/w{0), can be re- 
lated to the first moment f{do) in the limit ipi and 
ti2 oo. 

This formal treatment has been performed for a hop- 
ping model with discrete sites. This may be considered as 
an appropriate model also for more realistic systems for 
which the particles will fluctuate around the individual 
sites. In particular this is the case for the ion conduc- 
tor, studied in this work (see below for more details). In 
contrast, for glass-forming systems hopping dynamics is 
not so relevant. Therefore it is not possible to formulate 
such a simple relation between the derivative w{t) of the 
mean square displacement and the first moment f{do). 

Eq.|| also implies that the first moment f(do) in the 
limit toi does not contain new information as com- 
pared to the mean square displacement. For finite toi, 
however, the first moment can no longer be predicted 
from w(t). Thus new information as compared to the 
mean square displacement about the nature of correlated 
back and forth dynamics becomes available. Of partic- 
ular interest is the dependence of the first moment on 
rgi. For example one can learn whether there exists a 
long-range back-dragging effect, as implied e.g. in the 
percolation approach of ion dynamics, or whether after 
jumping to the next nearest neighbor site the memory 
about the initial site has been basically wiped out. 



IV. MODEL CALCULATIONS 

In order to clarify the information content of the first 
moment f{roi) and the second moment v{roi) we first 
calculate them for some simple one-dimensional models. 
We always consider the case of stochastic dynamics since 
we are interested in time scales which are beyond the 
ballistic regime. 

A. Harmonic oscillator 

This model is relevant to describe the backward and 
forward dynamics in the individual potential wells. We 
consider a harmonic oscillator with minimum at r = 0. 
Both times toi, ti2 are longer than the equilibration time 
of the oscillator, i.e. the mean square displacement is al- 
ready constant. We first calculate the probability pi(ri) 
that after the first time interval the system is at ri af- 
ter a motion of roi to the right. This can be formally 
calculated as 



Piin) 



dro po{ro)piri\ro)6{roi - (ri - rg)) (6) 
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Figure 2: f(roi) as obtained from Monte Carlo simulations of 
the simple periodic cos-potential. The potential is shown in 
the inset. 



Figure 3: f(roi) as obtained from Monte Carlo simulations of 
a periodic potential with alternating barriers. The potential 
is shown in the inset. 



where po('^o) denotes the equilibrium distribution and 
p(''iko) is the probability to move from rg to ri dur- 
ing time ^01- For large igi the latter term is identical 
to the equilibrium distribution, i.e. poiri). The result- 
ing gaussian integral can be easily solved. Here we are 
particularly interested in the first moment (ri) of pi(ri). 
One obtains 

in) = (l/2)roi. (7) 

For long ti2 the system acquires the average (r2) = 0, 
yielding 

r{roi) = {r2) - in) = -(l/2)roi. (8) 

Thus the back-dragging effect in the second time interval 
is proportional to the distance moved in the first time 
interval. 



ti2 — 10 tQi. The result of this simulation is shown in Fig. 
^. The dependence of ^(t'oi) on tqi can be understood 
from simple arguments: a) For rpi <C 1, one basically has 
the behavior of a harmonic oscillator, b) For tqi = 0.5, 
the particle has typically moved to a position close to the 
saddle between two wells. For infinite ioii simple symme- 
try considerations show that (ri) is exactly on the saddle. 
Since a particle on a saddle does not experience any ef- 
fective net force to any side one has f(0.5) = 0. For finite 
toi , a motion of rgi — 0.5 on average leaves the particle in 
the initial well |Q. Thus the function r(roi) has its zero 
for roi slightly larger than 0.5. c) For tqi approaching 
1, the particle has definitely crossed a saddle. Now the 
effective force points in the same direction as the initial 
jump direction. This results in a positive value of r(roi). 
d) For roi = 1, one has the same result as for roi = 0, 
i.e. f(roi) = 0. 



B. Periodic potential 

In the next step, we analyse a potential which involves 
local vibrations as well as hopping dynamics. Here we 
consider a periodic potential Epot{r) with minima at in- 
teger values of r, defined as 

F(r) = Fo[l - cos(27rr)]. (9) 

We are interested in the stochastic dynamics of a particle. 
In general, one has to resort to numerical simulations to 
calculate r(r'oi). Here we have modelled the dynamics via 
standard kinetic Monte Carlo simulations at the temper- 
ature T — 0.4Vo with step sizes much smaller than the 
distance of two adjacent minima. The qualitative fea- 
tures of the result do not depend on the exact value of 
this temperature. We have chosen toi such that on this 
time scale a particle leaves the initial well with a probabil- 
ity of approximately 50%. Furthermore we have chosen 



C. Potential with alternating barriers 

So far we have only discussed backwards dynamics due 
to intrawell dynamics. For an ion conductor one expects 
that dynamic forward backward correlations either result 
from static or from dynamic disorder. Here we briefly 
discuss a very simple model which contains non-trivial 
forward backward dynamics. It is shown in the inset of 
Fig. ||. Due to the alternating barrier heights a particle 
performs several forward and backward jumps until it can 
escape from the local cage. We denote the two different 
transition rates by Fi and r2 (Fi ^ F2). Our goal is to 
calculate r(roi — 1). For very long toi, i.e. toiT2 ^ 1 
the dynamics resembles that of a random walk so that 
forward backward correlations should not be relevant, i.e. 
f(l) w 0. Therefore we restrict ourselves to the case 
toil^2 ^ 1- The analysis is presented in Appendix |A|. 
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We obtain 
f(l) = - 

f(l) 



-(i/2)(i-2r2/ri)(i-2r2/ri) (lo) 
*(i-exp(-2riti2)) %ri<i 
-(i/2)(i-2r2/ri)(i-r2%) (ii) 
*(i-exp(-2riti2)) %ri>i. 



Thus there exists a Hmiting value for toi ~^ which is 
reached for ioi ~ 1/ri, i.e. the time scale of the fastest 
jump process. For larger values of toi the back-dragging 
effect decreases with time; see Eq.^ f(l) approaches 
zero for ^oi of the order of l/r2. The physical reason is 
that for longer times ipi the particle may also cross the 
high barrier during the first time interval. These events 
strongly reduce the total back-dragging effect in the sub- 
sequent time interval. Thus the toi-dependence contains 
valuable information about the time-scales involved in 
the dynamics and indicates at which time scale (here: 
^01 ~ 1/1^2) a simple random- walk description becomes 
relevant. 

For ti2 one obtains f — s- 0. This limit is trivial 
since there is no dynamics during the second time in- 
terval. In contrast, for ti2 00 the backjump effect is 
largest. Thus it is this limit which is relevant to judge 
the maximum backjump capabilities. The rest of the dis- 
cussion for this model system deals with this case. 

To show the full dependence of r(roi) on tqi we 
again performed Monte Carlo simulations for a poten- 
tial with alternating barriers. This potential was gen- 
erated from the cos-potential, discussed above, by scal- 
ing the cos-potential by a factor 1.5 in the intervals 
[-1,0], [1,2], [3,4], .... The temperature was T = 
0.51^0- ^01 is chosen such that the short-time limit 
toi < l/Fi is fulfilled whereas ti2 corresponds to the long- 
time limit ti2 ^ I/F2. The transition rate is propor- 
tional to the attempt frequency in the minimum, which 
scales with the square root of the force constant and 
with the Boltzmann factor. Thus one expects F2/F1 « 
N/L5exp((2 - 3)/0.5) « 0.17 which, according to EqjlO, 
yields f(l) « 0.22. 

The result for f(r'oi) is shown in Fig. |[ It resembles 
that of a periodic potential except for a systematic down- 
ward trend. Thus we have a superposition of the corre- 
lation effects of simple periodic potentials with wells and 
barriers and of barriers with different heights. Actually, 
it turns out that our estimate of f(l) agrees reasonably 
well with the numerical value of ca. 0.25. 



D. Random trap model 

Finally, we want to present a very simple model which 
allows one to grasp the relevant features of the sec- 
ond moment v{roi)- We start with an (possibly high- 
dimensional) array of traps with different depths. Such 
models have been extensively studied in the context of 
supercooled liquids |2^, ^ |3l|. Different rates F^ are 
randomly attributed to the different traps. For such a 
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Figure 4: V{t) as obtained from Monte Carlo simulations of 
a potential with random traps. The potential is shown in the 
inset. 



simple model the dynamics is purely diffusive. Here we 
are specifically interested in the ratio V{t) = w(roi — 
1, t)/v{roi = 0, t) for t^i = ti2 = t. According to our dis- 
cussion in Section III V{t) is a measure for the relevance 
of dynamic heterogeneities. 

As shown in Appendix ^ one can derive the relation 



y(i^O) = (F)(i 



(12) 



Without dynamic heterogeneities it does not matter 
whether or not a particle moves in the first time in- 
terval such that w(roi = l,i) = w(roi = 0,<). In case 
of no dynamic heterogeneities, i.e. a single value of F, 
one trivially has (F) = f/(f/F) and thus V(t) = 1. 
For a distribution of jump rates the product (F)(l/F) 
is larger than one. This can be easily rationalized for a 
bimodal rate distribution with rates Fi, F2 and weights 
fli, 02, respectively. A straightforward calculation yields 
V{t) = 1 + (0102) (Fi - F2)V(FiF2) which for a bimodal 
distribution is strictly larger than one. For stronger dy- 
namic heterogeneities, i.e. a broader distribution of rates 
F, V{t) also increases. Thus V{t ^ 0) is a direct measure 
for dynamic heterogeneities. 

In the opposite limit t — *■ 00, a particle with roi ~ 
has by chance returned to the original position. This im- 
plies that the condition roi = no longer implies any 
dynamical selection of slow particles. Thus one expects 
V{t 00) = 1. Whereas the detailed time-dependence 
of V{t) depends on more details of the model like the 
number of neighbor traps, the limiting values are gener- 
ally valid. 

In order to visualize the full time-dependence, and to 
check our analytic expression we have performed Monte- 
Carlo simulations for a one-dimensional random trap 
model. We have chosen two escape rates, characterized 
by ai = 0.035,02 = 0.965, Fi = 0.005, F2 = 0.1. For this 
specific choice of parameters one has (F)(l/F) — 1.61. 
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Figure 5: The derivative of the mean square displace- 
ment shown for different temperatures. Note that for t > Ips 
the function w{t) is mainly governed by long-range dynamics 
(see text). 
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Figure 6: The incoherent scattering function S{qmax,t) at 
different temperatures. The solid lines correspond to KWW- 
fits with P = 0.45. 



The time dependence of V{t) is shown in Fig. ^ as ob- 
tained via Monte Carlo simulations. For this simple 
model the algorithm can be implemented in a straightfor- 
ward way. One observes that the theoretical short-time 
and long-time limits are confirmed by the numerical data. 
Interestingly, V{t) displays a maximum. Qualitatively, 
this means that at the time scale of the maximum the ef- 
fects of dynamic heterogeneities are most pronounced. A 
more detailed discussion of the time-dependence and of 
the maximum is beyond the scope of the present paper. 



RESULTS 



Previous results 



It has been observed that lithium trajectories can be 
described as a series of local vibrations andjumps be- 
tween adjacent ionic sites; see e.g. Refs.[^ ^ for simi- 
lar features in previous simulations on alkali silicates. In 
our recent work we have shown that the van Hove self 
correlation function displays a strong peak for do « 2.6 
A. This peak is separated by a minimum at 1.5 A from 
the peak at the origin. The interpretation is straightfor- 
ward. The potential energy landscape, as supplied by the 
network, provides lithium sites with an average distance 
of 2.6 A which are separated by a saddle. 

Furthermore, it turned out that the mean square dis- 
placement curves at different temperatures show time- 
temperature superposition. In Fig. ^ we show their 
derivative, i.e. w[t)j for the five temperatures analysed 
in our prior work 0. As shown in the function w{t) 
for t > Ips is due to processes which involve long-range 
dynamical processes (|r(t) — r(0)| > 1.5 A). Whether or 
not these processes can always be interpreted as jumps is 
currently under investigation. For t < Ips the function 
is dominated by localized processes of the lithium 
ions. Since the presence of back and forth dynamics is 
equivalent to a decrease of w(t) with time, one directly 
sees that at least in the case of the two lowest tempera- 
tures T = TOOK and T = 640K long-range back and forth 
correlations are indeed important. In order to character- 
ize the typical time scales of the long-range dynamical 
processes we show in Fig. ^ the incoherent scattering 
function S{qmax,i) where qmax = 27r/do- It is a mea- 
sure of the probability that an ion is still (or again) at 
the initial site after time t. All data can be consistently 
fitted with a KWW-function f{t) oc ex.-p{-{t/TY) with 
13 ^ 0.45. 



B. Back and forth dynamics 

A quantitative analysis of the back and forth dynamics 
is possible on the basis of the first moment f(roi). In 
Fig. 0, we have plotted this function for a fixed value of 
ioi = 102.4 ps and different ti2 at T = 750 K. For large 
roi the statistics of these curves becomes quite poor since 
only a few ions participate. This problem becomes worse 
for short <oi and/or long ti2. In this and the following 
plot we restrict ourselves to the tqi -regions which possess 
a reasonable signal-to-noise ratio as estimated from the 
fiuctuations of the curves. 

On a qualitative level the dependence on Tqi resembles 
that shown in Fig. ^. For small roi one recovers the har- 
monic behavior, as seen from the very good agreement 
with the broken line r(roi) = — (l/2)roi. Furthermore 
the non-monotonic behavior directly reflects the presence 
of a saddle between adjacent lithium sites. Fig. |7| clearly 
reveals that the fraction of back and forth dynamics in- 
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Figure 7: The first moment r(roi) for toi = 102. 4ps and for 
different ti2 at T = 750 K. The broken line corresponds to 
f(roi) = -(l/2)oi. 
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Figure 8: The first moment r(do = 2.6A ) in dependence of 
toi ■ The value of ti2 has been scaled for every temperature on 
the basis of the diffusion constant [ti2(T = 640 K) = 3276.8 
ps]. In the inset also the values of toi are scaled by the diffu- 
sion constant such that ^ol(^ = 6407^) = t5i(r = 640is:). 



neighbor distance in the first time interval. We choose 
a large but constant value for ti2 (4.1 ns ps for T = 
640 K) and vary ioi- Comparison of different tempera- 
tures is achieved by choosing the respective value of ti2 
approximately proportional to the inverse diffusion con- 
stant. This results in tu = 26ps{T = 12A0K),ti2 = 
102ps{T = 980K),ti2 = 0.8ns{T = 750K),ti2 = 
2.0ns{T = 700K),ti2 = 4.1ns(T = 640K). The data 
are shown in Fig. |^. Obviously, for longer times toi the 
back-dragging effect becomes much smaller. This agrees 
with the theoretical considerations, discussed for the al- 
ternate barrier model. The limiting value f{do) = —1.3 A 
is only approximately reached for the lowest temperature 
T — 640 K. Interestingly, the dependence on toi is very 
gradual and extends over several decades of time. This 
shows that there exists a broad distribution of barriers 
experienced by the lithium ions. 

For a fixed value of toi the relevance of correlated back 
and forth dynamics increases with decreasing tempera- 
ture. From our previous discussion the short-time limit 
of toi « 1 ps is related to the dispersion according to Eq. 
^. Since only for toi > 1 ps long-range processes become 
relevant one should choose toi = 1 ps as the short-time 
limit in Eq||. Checking it, e.g. for T — 750K, one obtains 
f(do) = — 1 A which according to Eq. || corresponds to 
w{t — !■ oo)/w{lps) = 0.23. This agrees with the mea- 
sured value of ca. 0.2 (see Fig. ||). Note that Eq.^ is 
based on a strict hopping picture. Therefore one would 
not expect an exact agreement between both values. In 
any event, since only for T < 750 K f{do) comes close 
to the limiting value of —1.3 A the back-dragging ef- 
fect and, equivalently, the dispersion in the mean square 
displacement are relevant only in this low-temperature 
range. 

In the inset of Fig. || the individual curves are scaled 
with the scahng factor D{T)/D{T = 640X), thereby in- 
troducing the scaled time tgj^ . Within the statistical noise 
one can see a decent superposition, in particular for the 
three lowest temperatures. Thus one observes a time- 
temperature superposition also for this rather involved 
quantity of correlated back and forth dynamics. 



creases with increasing ti2. This is expected from our 
theoretical considerations; see above. 

Interestingly, for the two largest values of ti2 the func- 
tion r(roi) decays further for roi > do — 2.6 A. Thus 
back-dragging effects become stronger when jumping into 
the second nearest neighbor shell during toi. This obser- 
vation already goes beyond a scenario which is only based 
on back and forth correlations between adjacent lithium 
sites. 

Of particular interest is the dependence on tpi as al- 
ready discussed for the alternate barrier model. Since 
we are mainly interested in correlated back and forth dy- 
namics between nearest-neighbor positions we focus on 
the value of f(do) which characterizes the subsequent dy- 
namics of a particle which has jumped to the nearest 



C. Dynamic heterogeneities 

Finally we present results for the second moment 
i;(roi). In Fig. || we have plotted v{roi) for T — 750 
K and different values for toi = ti2. One can clearly dis- 
tinguish two regimes which are approximately separated 
by 2 A. For roi <ti 2 A , the second moment v{roi) is sig- 
nificantly smaller than for values tqi > 2 A. As discussed 
above this is a clear signature of dynamic heterogeneities. 
After a jump of 2 A an ion has basically achieved to cross 
the saddle and is already part of the adjacent ionic site. 
This means that a particle which has crossed the local 
saddle during the first time interval is much faster dur- 
ing the subsequent time interval. Starting from roi = do 
the second moment ^(roi) further increases with increas- 
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Figure 9: v{roi) atT— 750K for different choices of toi ~ ti2- Figure 10: V{t*) — v{dQ)/v{0) at different temperatures using 

scaled times. 



ing roi. This further increase, however, is much weaker. 
This observation would be compatible with only minor 
spatial correlations among ionic sites of similar mobility. 
In particular, it would contradict the scenario of compact 
regions of ionic sites each with a different ionic mobility. 
In this particle which has jumped twice would, on 

average, belong to a faster region than particles which 
only jumped once during the first time interval. This 
difference would show up in the mean square displace- 
ment during the second time interval since particles in 
the faster region would, on average, also jump further in 
the second time interval. This would lead to a strong 
increase of w(roi) beyond the nearest neighbor distance 
do. 

In order to study the temperature dependence and 
the time-dependence in greater detail we have calculated 
V{t) = v{do)/v{Q) in dependence of t toi ~ ti2 and for 
different temperatures. V{t) is a direct measure for the 
degree of dynamic heterogeneity on the length scale of the 
nearest-neighbor distance. In order to compare the dif- 
ferent temperatures we have again scaled all times by the 
ratio D{T)/D{T = 640K). The results are shown in Fig. 
[To| . One can clearly see that the degree of heterogeneity 
changes with temperature. The lower the temperature 



the stronger the dynamic heterogeneity. Thus on the 
level of dynamic heterogeneities the time-temperature 
superposition principle does not hold. For the lowest 
temperature V{t) is approximately 5.5 at the maximum. 
This number directly implies that the mean square dis- 
placement of a particle (corrected for possible backwards 
correlations) is 5.5 times larger in the second time inter- 
val if it has performed a jump in the first time interval 
as compared to particles which are still at the initial site 
after tgi- 



VI. DISCUSSION 

In this paper we have shown how analysis of three- 
time correlations can be used to get model-free informa- 
tion about the nature of the complex ion dynamics. The 
characteristics of back and forth dynamics is reflected by 
the first moment of the three-time conditional probabil- 
ity function, the dynamic heterogeneities by the second 
moment. The main results are: (i) the long-range back- 
ward correlations beyond the nearest neighbor position, 
(ii) the gradual decrease of backward correlations with 
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increasing toi, (iii) the time-temperature superposition 
principle for correlated back and forth dynamics, (iv) the 
significant dynamic heterogeneities at low temperatures, 
and (v) the lack of time-temperature superposition for 
the dynamic heterogeneities. In what follows we discuss 
these results in more detail. 

For ion conductors at low temperatures the basically 
immobile network serves as a pseudo-external field for the 
lithium ions. Therefore it has been attempted to model 
the lithium dynamics by simple lattice models like the 
random barrier or the random energy model [|[ |3^ . Also 
in these models the long-range backward correlations are 
present. Only in simple models like the alternate barrier 
model, discussed in this paper, back and forth jumps are 
restricted to the nearest neighbor position. The reason 
for these long-range backward correlations is quite intu- 
itive: the ions look for paths which can can be accessed 
rather easily. In particular at low temperatures these 
paths can be rather extended. Only for long times the 
ions manage to escape such a local path. This picture is 
consistent with the percolation approach which has been 
successfully applied to describe the dynamics in the ran- 
dom barrier and the random energy model 1^, |3^ . 

It may be interesting to compare this scenario with 
that of supercooled liquids. There all particles move 
on the same time scale. No pseudo-external fields are 
present. The complexity of the dynamics is due to the 
necessity of cooperative dynamics of the strongly inter- 
acting particles. A simple picture is to view a particle 
localized in the cage of the adjacent particles. The rel- 
evant relaxation process is to escape this local cage and 
afterwards being trapped in a new cage. In this case the 
first moment r{roi) is constant for values of rpi larger 
than the typical nearest-neighbor distance ||2l[. This is 
due to the fact that after leaving the initial cage no sig- 
nificant memory to that cage is left. The roi-dependence 
of r{rQi) for tqi > do shows that this simple cage picture 
cannot be used for the ionic dynamics. A possible ex- 
planation for this effect is the relevance of the static dis- 
ordered energy landscape in ion conductors which leads 
to a dramatic reduction of multi-particle correlations |^ 
such that cages, formed by adjacent ions, are less rele- 
vant. Alternatively, one might argue that the presence of 
long-range Coulomb interaction gives rise to long-range 
backward correlations. This aspect still has to be clari- 
fied. 

In the alternate barrier model the dependence of r{do) 
on ^01 refiects the values of the lowest and highest rele- 
vant barriers present in the system. On the time scale for 
which the highest relevant barrier is crossed no back and 
forth dynamics should be visible. Indeed, we see in Fig. 
^, e.g. for T = 750 K, that the back and forth dynamics 
becomes small for toi of the order of 1 ns. The relaxation 
time, i.e. the decay time of S{qmax, t), is of the same or- 
der. Within the alternate barrier model one expects that 
for tm somewhat smaller than l/r2 but still larger than 
l/Fi one expects a linear dependence of f(do) on t^i in 
marked contrast to the numerical results. This clearly 



shows that in contrast to the alternate barrier model the 
lithium metasilicate system is characterized by a broad 
distribution of relevant barriers. Furthermore no short- 
time limit is visible in Fig. || but f((io) decreases also 
for the shortest times which could be analysed. This 
indicates the presence of transitions with very small sad- 
dles or, possibly, broad anharmonic potentials which are 
already relevant on the time scale of a few hundred fem- 
toseconds. 

The alternate barrier model differs from our lithium 
metasilicate system also in another respect. For the alter- 
nate barrier model the typical jump time is of the order of 
1 /Fi. One can show that (i) the decay of w{t) is strongest 
for exactly this time scale and (ii) S{q — 2tt, t 1/Fi) is 
already significantly smaller than unity. For the lithium 
metasilicate system at, e.g., T — 640 K the decay of w{t) 
(due to long-range processes; see above) is strongest in 
the ps range whereas S{q,t) is still close to the short- 
time plateau value. This discrepancy could be alleviated 
by introducing an additional variation of site energies. 
This would shift the decay of S{q, t) to times of the order 
of I/F2. On a qualitative level this would imply that the 
short-time dynamics occurs in asymmetric double-well 
potentials. It remains an important question whether 
the same scenario also holds in lithium metasilicate, i.e. 
whether fast back and forth jumps in strongly asymmet- 
ric double- well potentials are present. Another scenario 
to rationalize the above-mentioned discrepancy will be 
presented further below. 

In Ref. 1^ it has been reported that for sodium silicate 
at T = 2000 K no backward correlations were observed. 
This is compatible with the experimental data for sodium 
silicate where the dispersion disappears for T « 1000 K. 
[5]. For the present system the dispersion disappears 
around a similar temperature. This shows up in a very 
little time-dependence of w{t) at 980 K in Fig. 5 beyond 
1 ps. Nevertheless, even these weak backward correla- 
tions are directly visible by analysis of the three-time 
correlations in Fig. 8. As discussed below these backward 
correlations may be due to a small number of ions. Thus 
a less detailed analysis may oversee these backward cor- 
relations. 

The approximate time-temperature superposition 
principle as seen in Fig. ^ implies that the back-dragging 
effect remains the same if analysed on appropriately ad- 
justed time scales. Note that this statement goes beyond 
the previous result that the mean square displacement 
displays time-temperature superposition. As discussed 
above the latter observable is only related to the ioi 
limit of r(roi). The time-tcmpcrature superposition in 
Fig. ^ may be used to discriminate between different 
models of ion dynamics. 

The significant heterogeneities, as characterized by the 
second moment, is compatible with the above-mentioned 
broad distribution of relevant barriers. Actually, such 
a distribution complicates the interpretation of the dis- 
persion w(t oo)/w{0) in terms of correlated back and 
forth dynamics. This is exemplified for the simple case of 
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two (temporarily distinct) ionic species with jump rates 
Fi ^ r2 which are present with probabihties pi and P2, 
respectively. If we have in mind a log-gauss distribution 
of rates one should choose pi > p2- Generalizing Eq.^ we 
get 



w{t ^ oo) 1 1 PiT -\- P2T 2r'^ 

w{t^{)) 2 rfo Pl^l+P2^2 



(13) 



The terms TiPi imply that during the short time interval 
toi the probability of a particle of species i to perform a 
hop is proportional to the rate and its occurrence proba- 
bility. In the limit, discussed above, and under the addi- 
tional assumption that |f^| is not much smaller or even 
larger than |f^|, this can be approximated by 



w{t 



00 



w(t 0) 



r 

da' 



(14) 



This result shows that the dispersion is to a large de- 
gree determined by the backjump properties of the fast 
species. Thus in case of significant dynamic hetero- 
geneities one has to be careful to relate the dispersion of 
the mean square displacement to the average backjump 
properties (here: {pif^ +P2r^)/do) of the ions. 

This observation allows one to envisage another ex- 
planation of the very different time scales where the de- 
cay of wit) is maximum (ps regime) and where the inco- 
herent scattering function S{qmax,t) decays (ns regime); 
see above. A few fast ions with significant back- and 
forth correlations can dominate the time-dependence of 
w{t) but hardly contribute to the decay of the incoher- 
ent scattering function. Actually, for the random energy 
model we have observed that the back and forth corre- 
lations are strongest for the fast particles such that this 
effect is indeed present there. Whether it is this sce- 
nario or the above-discussed presence of energetic disor- 
der, which dominates the different time dependence of 
w{t) and S{qmax,t) for lithium metasilicate, awaits fur- 
ther clarification. 

We have shown that the degree of heterogeneity de- 
pends on temperature. In contrast, the incoherent scat- 
tering function S{qmaxTt) fulfills the time-temperature 
superposition principle. Actually, the same observations 
have been made for glass- forming liquids . For these 
systems the non-exponentiality of S{qmax,t) mainly re- 
flects the broad distribution of relaxation times. Thus 
naively one would expect a lower value of the KWW- 
exponent /3 with increasing degree of heterogeneity. It 
still has to be shown why at least in the temperature 
range, accessible to simulations, this is not the case. 

It may be interesting to compare our approach with 
that of Habasaki and Hiwatari For the same sys- 

tem at T = 700 K they determined during 1 ns the dis- 
tribution of square displacements of all particles. They 
observed a large variance, thus indicating some distribu- 
tion of relaxation times. On a qualitative level this result 
is compatible with the results reported above. One ad- 
vantage of the present approach is that via study of the 



roi-dependence also information about length scales are 
available. Furthermore, the value of V{t — > 0) has a di- 
rect interpretation in terms of local rate distributions. 

Having identified several properties of the complex ion 
dynamics in quantitative terms one would like to relate 
them to more microscopic properties like the distribution 
of oxygens around the lithium ions and to see whether 
mainly the interaction among the different ions or the 
interaction with the basically static network gives rise to 
the observations reported in this work. 
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Appendix A: 

Calculation of r(roi = 1) proceeds in two steps. First 
we calculate the probability that a particle has moved 
by one unit to the right during some time ioi- Since we 
have two non-equivalent sites we have to calculate, on 
the one hand, the transition probability q^o (eo stands 
for even-odd) from r = to r = 1 and, on the other 
hand, the transition probability qoe from r = 1 to r = 2. 
In the second step we calculate the average motion Ar^ 
(starting, e.g., from r = 2) and Avg (starting, e.g., from 
r — 1) during the second time interval. One expects 
Are > and Atq < and for reasons of symmetry Ar^ — 
—Aro which we abbreviate as Ar. With this information 
we can finally calculate 



r(l) 



Aroqeo + Areq, 



oe * qeo 

= —Ar 



9o 



(Al) 



For the calculation of qeo and qoe we take into account 
that in the limit toir2 <C 1 multiple transitions over the 
higher barrier can be neglected. We start with a particle 
either at r = or r = 1 and we are interested in the prob- 
ability to be at r = 1 or r = 2 after time toi 1 respectively. 
Neglecting those terms, which are only relevant in case 
that a particle has crossed a high barrier at least twice, 
we end up with the following system of rate equations for 
the site populations pi 

{d/dt)po = -(ri + r2)po + ripi (A2) 
id/dt)pi = -(ri + r2)pi + ripo (A3) 
id/dt)p2 = -riP2 + riP3 + r2Pi (A4) 

id/dt)p3 = -Tip3 + Tip2. (A5) 

This set of differential equations can be directly solved. 
qeo can be identified as pi with initial condition r = 
and qoe as p2 with initial condition r = 1. One obtains 
after a short calculation 

qeo = (1/2)(1 - r2toi)(l - exp(-2Fitoi)) (A6) 
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and 



qoe = 7fT-[l - exp(-2ritoi) + ritoi(l + exp(-2ritoi)]. 
41 1 



Thus we finally get (using again Fi 

^!^-Z«££ = 1 - 2^^ (toiF « 1) 

9eo + 9oe 1 1 



(A7) 



(A8) 



and 



9. 



eo Hoe 



Qeo + qoe 



= 1 - Tatoi (toiT > 1). (A9) 



As a second ingredient we want to calculate Ar which 
is the average coordinate (r) after time ti2 for the initial 
condition r = 0. Using the standard trick of introducing 
the functions 



and 



Tq = ^{-lypjexpiigj) 



(AlO) 



(All) 



one can write down two linear differential equations in- 
volving Tq and Sq 



which in the limit Fi 3> T2 can be rewritten as Ar = 

(i/2)(i - 2r2/ri)(i - exp(-2riti2)). 

Thus we have calculated all ingredients which are nec- 
essary for determination of f(l). 



Appendix B: 

Here we calculate the short-time limit of the random 
trap model. We always choose toi — ti2 = t. Wc define 
the probability that a trap has the escape rate Fj as Pi. 
In equilibrium the probability that a particle is in a trap 
with escape rate Fj is proportional to Pi/Fj. Then we 
can write 



v{roi = l,t) 



(Bl) 



Here qi{t) denotes the probability that after time t the 
particle has moved by one unit, starting in a trap with 
escape rate Fj and r| {t) the short-time expression for the 
mean square displacement, starting from a trap with rate 
Tj. One simply has qi(t) ~ Tit and r'j{t) = Tjt. After a 
short calculation one obtains 

v{roi =l,t^O)= Y,Pi^i = (r)- (B2) 



{d/dt)Sq{f) = -(Fi + T2)Sq{l - cosq) - i(Fi - F2)7^Ailfi^ 
{d/dt)Tq{t) = -(Fi+F2)Tq(l+cosg) + «(Fi-F2)5^MB^ 

which can be solved with standard methods after spec- 
ification of the initial condition. Of interest for us is the 
expectation value (r) for the initial condition r(0) = 0. 
This expression can be calculated from 



In analogy one obtains for the short-time expansion of 



(r) = —i \im{d/dq)Sq. 

q^O 



(A14) 



It turns out 



Ar = (l/2)^i-^ (1 - exp(-2(Fi + F2)ii2)) (A15) 
J- 1 + i 2 



v{roi = 0, t) 



Thus one obtains 

v{roi = l,t) 
v{roi = 0,t) 



(B3) 



(F)(l/F). 



(B4) 
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